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ABSTRACT 

Prospects for SO(10) as a minimal grand unification group have recently 
been heightened by several considerations such as the MSW resolution of the so- 
lar neutrino puzzle, baryogenesis, possibility for understanding fermion masses etc. 
I review the present status of the minimal SO(10) models with special emphasis 
on the predictions for proton lifetime and predictions for neutrino masses for the 
non-supersymmetric case and discuss some preliminary results for the supersym- 
metric case. It was generally believed that minimal SO(10) models predict wrong 
mass relations between the charged fermions of the first and second generations; 
furthermore, while the smallness of the neutrino masses in these models arises from 
the see-saw mechanism, it used to be thought that detailed predictions for neutrino 
masses and mixings require further adhoc assumptions. In this talk, I report some 
recent work with K.S.Babu, where we discovered that the minimal SO(10) model, 
both with and without SUSY, has in it a built-in mechanism that not only corrects 
the bad mass relations between the charged fermions but at the same time allows a 
complete prediction for the masses and mixings in the neutrino sector. We define 
our minimal model as the one that consists of the smallest set of Higgs multiplets 
that are needed for gauge symmetry breaking. Our result is based on the hypothe- 
sis that the complex 10 of Higgs bosons has only a single coupling to the fermions. 
This hypothesis is guaranteed in supersymmetric models and in non-SUSY models 
that obey a softly broken Peccei-Quinn symmetry. 

1. Introduction: 

At the present time ,the only experiments where there is some hint of new 
physics beyond the standard model are the ones which have detected the neutrinos 
emitted from the solar core[l]. The deficit of solar neutrinos reported in these 
experiments from Homestake, Kamiokande, SAGE, and GALLEX[2] if confirmed in 
the other proposed experiments such as the BOREXINO[3] and SNO[4] will confirm 
that the neutrinos have masses and mixings very much like the quarks. The observed 
deficit in the present experiments can be explained in terms of neutrino oscillations 
in two different ways: (i) long wave length vacuum oscillation [5], and (ii) resonant 
matter oscillation (the Mikheyev- Smirnov-Wolfenstein (MSW) effect [6]). Assuming 
a two-flavor u e - oscillation, in the former case, the neutrino masses and mixing 
angle should satisfy Am 2 ~ 10~ 10 eV 2 and sin 2 26> e/l ~ (0.75 to 1). In case of MSW there 
are two allowed windows that fit all of the experimental data[7]: (a) the small mixing 

1 Invited talk presented at the workshop on "New Physics with New Experiments" held in Kaimierz, 
Poland, May, 1993. 



angle non-adiabatic solution, which requires Am 2 ~ (0.3 to 1.2) x 10~ 5 eV 2 and sin 2 26> eAI ~ 
(0.4 to 1.5) x 10~ 2 , and (b) the large angle solution with Am 2 ~ (0.3 to 5) x 10~ 5 eV 2 and 
sin 2 26> eM ~ (0.5 to 0.9). In all these cases, barring an unlikely scenario of near mass 
degeneracy among neutrinos, either or v T should have mass in the (10~ 5 to 10~ 3 ) eV 
range. Specifically for the MSW resolution, the mass has to be of order 10~ 3 eF. 

There are two other indications, though more controversial, for a non- vanishing 
neutrino mass; i) the apparent deficit of GeV muon neutrinos in the cosmic rays that 
originate from the decay of kaons and pions[8] in the atmosphere; and ii) the need 
for a hot component in the dark matter of the universe[9]. The first effect would be 
an indication of a possible oscillation between the and v T , again indicating the 
existence of a non- vanishing mass for the neutrinos. The simplest candidate for the 
hot component of the dark matter is a neutrino wth mass in the few electron volt 
range. The controversy surrounding the first effect has to do with the uncertainties 
in the atmospheric muon neutrino fluxes [10], whereas that with the second is due 
to possible effects of a cosmological constant and other cosmological inputs that go 
into the discussion of structure formation. In fact if all these three effects are to be 
taken together, the resulting form of the neutrino mass matrix becomes very seri- 
ously constrained [11] and one has to invoke specific kinds of family symmetries to 
understand them. Our goal in this article will be to see to what extent the present 
observations motivate the serious consideration of SO(10) as a grand unification 
symmetry and how much of the data can be understood in this framework. 

A natural explanation for the origin of such tiny neutrino masses is the see- 
saw mechanism[12], which is based on a mass matrix of the following form: 



Here ,the m D and M N are 3x3 matrices denoting the Dirac and the Majorana 
masses involving the left and the right-handed neutrinos. The diagonalization of 
the see-saw matrix leads to the light neutrino masses of the following form: 



The mass matrix M N corresponds to the scale of B-L breaking and represents 
physics beyond the standard model. It therefore could be large whereas m D is 
characteristic of the electroweak scale and is in the GeV range, thus explaining the 
smallness of the neutrino masses. The solar neutrino puzzle indicates that the B — L 
scale is in the (10 12 - 10 16 ) GeV range. 

All of the observations above, viz., non-zero neutrino masses, the see-saw 
mechanism, and a high B-L scale, fit rather naturally in grand unified models 
based on the gauge group 5O(10). In its non-supersymmetric version, experimental 
constraints from proton life-time and the weak mixing angle sin 2 6 w require that 
5O(10) breaks not directly into the standard model, but at least in two steps. In 
a two-step breaking scheme, the left-right symmetric intermediate scale is around 
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10 12 GeV[13]. In supersymmetric 50(10) there is no need for an intermediate scale, 
50(10) can break directly to the standard model at around 10 16 GeV. 

There is also another compelling cosmological reason for the 50(10) model 
which is absent in GUT models , which has to do with the generation of cosmo- 
logical baryon asymmetry. It has been pointed out that, in 50(10) models with 
the see-saw mechanism ,the right-handed neutrino is very heavy and has a mass 
near the B - L breaking scale of 10 12 GeV or so. Therefore,its decay can generate 
a lepton asymmetry when the temperature of the Universe is about 10 11 GeV or 
so. This lepton asymmetry gets subsequently converted to the baryon asymmetry 
due the sphaleron effects at around T ~ 300GeV[14]. In investigating the survival of 
the lepton asymmetry from this high temperature down to the electro-weak phase 
transition temperature , one must make sure that the lepton violating interactions 
are out of equilibrium. This is guaranteed only if the B — L symmetry scale is bigger 
than 10 11 or so. This therefore again points towards an 50(10) theory where such 
high scales are naturally generated. 

To confront 50(10) models with the solar neutrino data, one must make pre- 
cise predictions of the neutrino masses and mixing angles. This requires, however, 
detailed information of the Dirac neutrino mass matrix as well as the Majorana 
matrix. In grand unified theories (GUTs), it is possible to relate the quark masses 
with the lepton masses. It used to be thought that in simple 50(10) models, the 
charge -1/3 quark mass matrix is related to the charged lepton matrix and the neu- 
trino Dirac mass matrix is related to the charge 2/3 quark matrix at the unification 
scale. No simple simple way was known, in general, to relate the heavy Majorana 
matrix to the charged fermion observables. This prevented any prediction of light 
neutrino spectrum without making extra symmetry assumptions. 

In a recent paper, [15] it was shown Babu and this author that this situation 
resulted from an incomplete analysis of the Higgs sector of the model and we showed 
that in a class of minimal 50(10) models, not only the Dirac neutrino matrix, but 
the Majorana matrix also gets related to observables in the charged fermion sector. 
This leads to a complete prediction for the neutrino masses and mixings without 
adding any extra Higgs multiplets or any new symmetries to the theory. We use a 
simple Higgs system with one (complex) 10 and one 126 that have Yukawa couplings 
to fermions. The 10 is needed for quark and lepton masses, the 126 is needed for 
the see-saw mechanism. Crucial to the predictivity of the neutrino spectrum is 
the observation that the standard model doublet contained in the 126 receives an 
induced vacuum expectation value (vev) at tree-level. In its absence, one would 
have the asymptotic mass relations m b = m T , m s = m^, m d = m e . While the first 
relation would lead to a successful prediction of m b at low energies, the last two are 
in disagreement with observations. The induced vev of the standard doublet of 126 
corrects these bad relations and at the same time also relates the Majorana neutrino 
mass matrix to observables in the charged fermion sector, leading to a predictive 
neutrino spectrum. 
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2. Minimal SO(10) GUT without supersymmetry: 



In this section, we shall consider non-SUSY 50(10) model to illustrate our 
mechanism. In the minimal scenerio defined here, the 50(10) symmetry breaks to 
the standard model via the SU{2) L x SU{2) R x SU(4) C = G 22 4 chain. The breaking of 
50(10) via G 2 24 is achieved by either a 54 or a 210 of Higgs. The 210 also breaks the 
discrete L>-parity,[16] whereas the 54 preserves it. .D-parity is a local discrete Z 2 - 
subgroup of 5O(10), under D, a fermion field / transforms into its charge conjugate f c . 
Breaking of D-parity at the GUT scale makes the see-saw mechanism natural [17] .It 
further changes the evolution of the gauge coupling constants by making the Higgs 
spectrum left-right asymmetric. It also eliminates the cosmological domain wall 
problem that can arise if this Z 2 local symmetry survived to the intermediate scale. 
We will therefore work with the 210 Higgs multiplet, although our mechanism to 
cure the fermion mass problem and predict the neutrino masses is independent of 
this choice. We will denote this chain with SU(A) C as intermediate symmetry as chain 
A. If instead of a 210 dim. Higgs multiplet, we chose a 45 plus 54 combination, then 
the 5O(10) symmetry would break down to SU{2) L x SU(2) R x U(1) B -l x SU(3) c = G 22 i 3 
group. In this case, also D-parity breaks at the GUT scale , thereby making the 
see-saw formula natural. We will denote this chain B. The second stage 

of symmetry breaking goes via the 126 in both cases. Finally, the standard model 
electro-weak symmetry breaking proceeds via the 10. 

2a. Mass scales and Proton decay: 

Before discussing the fermion sector of the theory, let us discuss the mass 
scales as predicted by the low energy values of sin 2 9 w , a strong and a em . This analysis 
was first carried out for chains with D-parity breaking by Chang et. al. in ref.13 
and recently been reanalyzed by Deshpande et al.[13]. The threshold corrections 
due to heavy particles was carried out recently by Parida and this author [18]. In 
the threshold correction analysis, the survival hypothesis was used to determine 
which Higgs submultiplets are at what mass scale and then allowing for reasonable 
uncertainties in the heavy Higgs boson masses, the uncertainties in the mass scales 
was estimated. It was found that the different mass scales and the uncertainties in 
their values depend on the symmetry breaking chain and we have: 



Chain A: 



M„ = 10 15 - 8 -i 8 ? ± - 2 GeV 



Mj = 10 



,11.5 



+ 2.8 
'-1.5 



±.02 



GeV 



(3) 



Chain B: 



M v = 10 



i 15.8±.l±.25 



GeV 



Mi = 10 9 - 3 ±18 Get/ 



(4) 
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The significance of these results is that the M v governs the life-time of the 
proton in 50(10) models whereas the Mj governs the neutrino masses vis the see-saw 
mechanism. We will come to the neutrino masses later. The proton life-time in 
both symmetry breaking chains is given by [18]: 

ChainA : t p = 1.6 x 10 35± - 7± - 9 - 3 ^ years (5) 

chainB : r p = 1.6 x 10 35± - 7±1 ' 0± - 8 years (6) 

In the above formulae, the last uncertainties are from the threshold cor- 
rections whereas the the first and second are from the matrix element and a str 
uncertainties respectively. It is clear from eq.(3) that the threshold uncertainties 
in Mi and M v are so large that one could consider this as an almost single scale 
theory. This possibility has been studied in detail in ref.19 and indeed, the results 
turn out to be quite consitent with the present lower limits on proton lifetime. 

2b. Fermion masses: 

Let us now turn to the fermion sector, we denote the three families belonging 
to 16-dimensional spinor representation of 50(10) by ip a , a = 1 - 3, the complex 10- 
plet of Higgs by H, and the 126-plet of Higgs by A, the Yukawa couplings can be 
written down as 

L Y = K b 4> a 4> b H + fabipaipbA + B.C. (7) 

Note that since the 10-plet is complex, one other coupling ip a ipb~H is allowed in 
general. In SUSY-5O(10), the requirement of supersymmetry prevents such a term. 
In the non-SUSY case, which we are focussing on in this section, we forbid this 
term by imposing a U(l) PQ symmetry, which may anyway be needed in order to 
solve the strong CP problem. 

The 10 and 126 of Higgs have the following decomposition under G 2 24'- 126 — > 
(1,1,6) + (1,3,10) + (3,1, TO) + (2,2,15), 10 -> (1,1,6) + (2,2,1). Denote the (1,3,10) and 
(2, 2, 15) components of A(126) by A R and S respectively and the (2, 2, 1) component 
of H(10) by $. The vev < A^ >= v R - 10 12 GeV breaks the intermediate symmetry 
down to the standard model and generates Majorana neutrino masses given by fv R . 
<f> contains two standard model doublets which acquire vev's denoted by n u and n d 
with K„ ;d ~ 10 2 GeV. k u generates charge 2/3 quark as well as Dirac neutrino masses, 
while n d gives rise to -1/3 quark and charged lepton masses. 

Within this minimal picture, if k u , n d and v R are the only vev's contributing 
to fermion masses, in addition to the SU(5) relations m b = m T , m s = m^, m d = ™ e , 
eq. (7) will also lead to the unacceptable relations m u : m c : m t = m d : m s : m b . 
Moreover, the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix will be identity. 
Fortunately, within this minimal scheme, we have found new contributions to the 
fermion mass matrices which are of the right order of magnitude to correct these 
bad relations. To see this, note that the scalar potential contains, among other 
terms, a crucial term 

Vi = AAAAff + H.C. (8) 
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Such a term is invariant under the U(1)pq symmetry. It will be present in the SUSY 
50(10) as well, arising from the 210 F-term or perhaps as a Planck induced term 
etc. This term induces vev's for the standard doublets contained in the S multiplet 
of 126. The vev arises through a term ~K R A R i:$ contained in V x . 

We can estimate the magnitudes of the induced vev's of £ (denoted by v u 
and v d along the up and down directions) assuming the survival hypothesis to hold: 

v u ,d ~ A (^T—j ^ ■ (9) 

Suppose M v ~ 10 15 GeV, Mi ~ 3 x 10 12 GeV as in the chain A and M s - 10 14 GeV, 
consistent with survival hypothesis, then v u and v d are of order 100 MeV, in the right 
range for correcting the bad mass relations. We emphasize that there is no need for 
a second fine-tuning to generate such induced vev's. In the SUSY version with no 
intermediate scale, the factor (v R /M£) is not a suppression, so the induced vev's can 
be as large as n Uid . This is the key observation of ref.15, which leads to the predictive 
5O(10) model in the neutrino sector, while removing the disagreement with charged 
fermion, masses that were thought to exist in the minimal model before. 

Another way to view this is to realize that the effect of the mixing is to leave 
a pair of light standard model doublets are admixtures of doublets in (2, 2, 1) and 
(2, 2, 15) multiplets and these light doublets acquire vevs at low energies leading to 
the same mass patterns as before. 

It is also worth noting that if we considered the symmetry breaking chain 
B, then since the intermediate scale is at most 10 10 GeV, the induced vevs v u and v d 
would be very small and would not provide any useful correction to the bad charged 
fermion mass relations. Thus, our scenerio would prefer the SU(4) C intermediate 
symmetry chain. 

3. The charged fermion masses: 

We are now in a position to write down the quark and lepton mass matrices 
of the model: 

M u = hn u + fv u M d = hn d + fv d 

M? = hn u - 3fv u M, = h Kd - 3fv d 

Mff = fv R . (10) 

Here Mf? is the Dirac neutrino matrix and M^f is the Majorana mass matrix. 

Before proceeding, we should specify the origin of CP violation in the model. 
We shall assume that it is spontaneous or soft, that will keep the number of parame- 
ters at a minimum. The Higgs sector described above already has enough structure 
to generate realistic CP violation either softly or spontaneously. The Yukawa cou- 
pling matrices h and / in this case are real and symmetric. Although there will be 
three different phases in the vev's (one common phase for k u and n d and one each 
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for v u and v d ), only two combinations enter into the mass matrices, as the overall 
phase can be removed from each sector. We shall bring these two phases into v u 
and v d and hence forth denote them by v u e %a and v d e il3 . 

To see the predictive power of the model as regards the neutrino spectrum, 
note that we can choose a basis where one of the coupling matrices, say h, is real 
and diagonal. Then there are 13 parameters in all, not counting the superheavy 
scale v R : 3 diagonal elements of the matrix hn u , 6 elements of fv u , 2 ratios of vev's 
r\ = Kd/n u and r 2 = v d /v u , and the two phases a and (3. These 13 parameters are 
related to the 13 observables in the charged fermion sector, viz., 9 fermion masses, 
3 quark mixing angles and one CP violating phase. The light neutrino mass matrix 
will then be completely specified in terms of other physical observables and the 
overall scale v R . That would lead to 8 predictions in the lepton sector: 3 leptonic 
mixing angles, 2 neutrino mass ratios and 3 leptonic CP violating phases. 

The relations of eq. (10) hold at the intermediate scale Mi where quark- 
lepton symmetry and left-right symmetry are intact. There are calculable renor- 
malization corrections to these relations below Mi. The quark and charged lepton 
masses as well as the CKM matrix elements run between Mi and low energies. 
The neutrino masses and mixing angles, however, do not run below Mj, since the 
right-handed neutrinos have masses of order Mi and decouple below that scale. The 
predictions in the neutrino sector should then be arrived at by first extrapolating 
the charged fermion observables to Mi. 

We fix the intermediate scale at M d = 10 12 GeV and use the one-loop standard 
model renormalization group equations to track the running of the gauge couplings 
between M z and Mj. 

To compute the renormalization factors, we choose as low energy inputs the 
gauge couplings at M z to be ai(M z ) = 0.01688, a 2 (M z ) = 0.03322, a 3 (M z ) = 0.11. For 
the light quark (running) masses, we choose values listed in Ref. 20. The top-quark 
mass will be allowed to vary between 100 and 200 GeV. Between 1 GeV and M z , we 
use two-loop QCD renormalization group equations for the running of the quark 
masses and the SU(3)c gauge coupling, [20] treating particle thresholds as step func- 
tions. From M z to Mi, the running factors are computed semi-analytically both for 
the fermion masses and for the CKM angles by using the one-loop renormalization 
group equations for the Yukawa couplings and keeping the heavy top-quark contri- 
bution[21]. The running factors, defined as % = m^M^/m^rrii) [ry, = m^M^/m^l GeV) 
for light quarks {u,d,s)] are r](u,c,t) = (0.273,0.286,0.506), r](d,s,b) = (0.279,0.279,0.327), 
r](e,fj,,T) = 0.960 for the case of m t = 150 GeV. The (common) running factors for 
the CKM angles (we follow the parameterization advocated by the Particle Data 
Group) S 2 3 and S 13 is 1.081 for m t = 150 GeV. The Cabibbo angle S 12 and the KM 
phase 5km are essentially unaltered. 

Let us first analyze the mass matrices of eq. (10) in the limit of CP con- 
servation. We shall treat spontaneous CP violation arising through the phases of 
the vev's v u e ia and v d e if3 as small perturbations. This procedure will be justified 
a posteriori. In fact, we find that realistic fermion masses, in particular the first 
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family masses, require these phases to be small. 

We can rewrite the mass matrices M U M® and M^f of eq. (10) in terms of the 
quark mass matrices and three ratios of vev's - n = n d /n u , r 2 = v d /v u ,R = v u /v R : 

r 2 -n r 2 - n 

M? = 3 -^M u A —M d , 

r 2 - n r 2 - n 

= ^^^M u -^—^—M d . (11) 
R n — r 2 R n — r 2 

It is convenient to go to a basis where M u is diagonal. In that basis, M d is 
given by M d = V M f agonal V T ', where M d d ia9 ° nal = diagonal(™ d , m„ m h ) and V is the CKM 
matrix. One sees that Mi of eq. (11) contains only physical observables from the 
quark sector and two parameters n and r 2 . In the CP-conserving limit then, the 
three eigen-values of Mi will lead to one mass prediction for the charged fermions. 
To see this prediction, Mi needs to be diagonalized. Note first that by taking the 
Trace of Mi of eq. (11), one obtains a relation for n in terms of r 2 and the charged 
fermion masses. This is approximately 

r\ ~ (m T + 3mb) /4m* (12) 

(as long as r 2 is larger than m b /m t ). Since \m b \ ~ |m T | at the intermediate scale to 
within 30% or so, depending on the relative sign of m h and m T , n will be close to 
either m b /m t or to (m h /2m t ). Note also that if r 2 > n, M ; becomes independent of r 2 , 
while M^ 5 retains some dependence: 

Mi ~ 4nM„ - 3M d , Mjf ~ Af u - —M d . (13) 

This means that the parameter r 2 will only be loosely constrained from the charged 
fermion sector. 

We do the fitting as follows. For a fixed value of r 2 , we determine n from 
the Tr(M/) using the input values of the masses and the renormalization factors 
discussed above. M t is then diagonalized numerically. There will be two mass 
relations among charged fermions. Since the charged lepton masses are precisely 
known at low energies, we invert these relations to predict the d-quark and s-quark 
masses. The s-quark mass is sensitive to the muon mass, the d-mass is related 
to the electron mass. This procedure is repeated for other values of r 2 . For each 
choice, the light neutrino masses and the leptonic CKM matrix elements are then 
computed using the see-saw formula. 

4. Predictions for neutrino masses and mixings: 

We find that there are essentially three different solutions. A two-fold ambi- 
guity arises from the unknown relative sign of m b and m r at Mj. Although solutions 
exist for both signs, we have found that a relative minus sign tends to result in 
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somewhat large value of m s /m d . Our numerical fit shows that the loosely con- 
strained parameter r 2 cannot be smaller than 0.1 or so, otherwise the d-quark mass 
comes out too small. Now, the light neutrino spectrum is sensitive to r 2 only when 
r 2 ~ Am s /m c - ±0.4, since the two terms in become comparable (for the sec- 
ond family) then. Two qualitatively different solutions are obtained depending on 
whether r 2 is near ±0.4 or not. 

Numerical results for the three different cases are presented below. The input 
values of the CKM mixing angles are chosen for all cases to be 5 i2 = -0.22, 5" 23 = 
0.052, 5i 3 = 6.24 x 10~ 3 . Since 8 KM has been set to zero for now, we have allowed 
for the mixing angles to have either sign. Not all signs result in acceptable quark 
masses though. Similarly, the fermion masses can have either sign, but these are 
also restricted. The most stringent constraint comes from the d-quark mass, which 
has a tendency to come out too small. Acceptable solutions are obtained when 
6*23, 0i3 are in the first quadrant and i2 in the fourth quadrant. 

Solution 1: 



Input : m u (l GeV) = 3 MeV, m c (m c ) = 1.22 GeV, m t = 150 GeV 
m b (m b ) = -4.35 Ge^ n = -1/51.2, r 2 = 2.0 
Output : m d (l GeV) = 6.5 MeV, m s (lGeV) = 146 MeV 
(m„ e , m„ M , m„ T ) = R (2.0 x 10~ 2 , 9.9, -2.3 x 10 4 ) GeV 

/ 0.9488 0.3157 0.0136 \ 



^lcpton 
KM 




(14) 



Solution 2: 



Input : m u (l GeV) 



3 MeV, m c {m c ) = 1.22 GeV, m t = 150 GeV 



m b (m b ) 



-4.35 Gey, n = -l/51, r 2 = 0.2 



Output : m d (l GeV) 
(m Vti ,m Vli ,m VT ) 



5.6 MeV, m s (l GeV) = 156 MeV 
i?(7.5 x 10" 3 , 2.0, -2.8 x 10 3 ) GeV 



^lcpton 
KM 




(15) 



Solution 3: 



Input :m„(l GeV) = 3 MeV, m c (m c ) = 1.27 GeV, m t = 150 GeV 



m b (m b ) = -4.35 Gey n = -1/51.1, r 2 = 0.4 



Output : m d {\ GeV) = 6.1 MeV, m s (l GeV) = 150 MeV 
{m Ve ,m^,m VT ) = R (4.7 x 10~ 2 , 1.4, -5.0 x 10 3 ) GeV 
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/ 0.9966 0.0627 -0.0541 \ 
V KM° n = I -0-0534 0.9858 0.1589 . (16) 
V 0.0633 -0.1555 0.9858 / 

Solution 1 corresponds to choosing n ~ m b /m t . All the charged lepton masses 
are negative in this case. Since r 2 is large, the Dirac neutrino matrix is essentially 
M u , which is diagonal; so is the Majorana matrix. All the leptonic mixing angles 
arise from the charged lepton sector. Note that the predictions for m d and m s are 
within the range quoted in Ref. 20. The ratio m s /m d = 22 is within the allowed 
range from chiral perturbation theory estimates [22]. The mixing angle sin^ c _^ 
relevant for solar neutrinos is 0.30, close to the Cabibbo angle. Such a value may 
already be excluded by a combination of all solar neutrino data taken at the 90% CL 
(but not at the 95% CL)[7]. Actually, within the model, there is a more stringent 
constraint. Note that the v^—v T mixing angle is large, it is approximately 3\V C b\ ~ 0.16. 
For that large a mixing, constraints from - v T oscillation experiments imply [23] 
that \ml T - ml \ < 4 eV 2 . Solution 1 also has m VT /m v ^ ~ 2.3 x 10 3 , requiring that 
< 0.9 x 10~ 3 eV. This is a factor of 2 too small for v e - MSW oscillation for 
the solar puzzle (at the 90% CL), but perhaps is not excluded completely, once 
astrophysical uncertainties are folded in. If v T mass is around 2 x 10~ 3 eV, v e - v T 
oscillation may be relevant, that mixing angle is ~ 3\V td \ — 6%. It would require the 
parameter R = v u /v R ~ 10~ 16 or v R - 10 16 GeV for v u ~ 1 GeV. Such a scenario fits 
very well within SUSY-SO(IO). 

Solution 2 differs from 1 in that r 2 is smaller, r 2 = 0.2. The ratio m s /m d = 27.8 
is slightly above the limit in Ref. 22. The 1-2 mixing in the neutrino sector is large 
in this case, so it can cancel the Cabibbo like mixing arising from the charged lepton 
sector. As we vary r 2 from around 0.2 to 0.6, this cancellation becomes stronger, 
the v e - mixing angle becoming zero for a critical value of r 2 . For larger r 2 , the 
solution will approach Solution 1. The - v T mixing angle is still near 3\V c b\, so 
as before, to„ t < 2 eV. From the v T /v^ mass ratio, which is 1.4 x 10 3 in this case, we 
see that < 1.5 x 10~ 3 eV. This is just within the allowed range[7] (at 95% CL) 
for small angle non-adiabatic v e - MSW oscillation, with a predicted count rate 
of about 50 SNU for the Gallium experiment. Note that there is a lower limit of 
about 1 eV for the v r mass in this case. Forthcoming experiments should then be 
able to observe - v T oscillations. A v T mass in the (1 to 2) eV range can also be 
cosmologically significant, it can be at least part of the hot dark matter. In SUSY 
50(10), v e - v T oscillation (the relevant mixing is about 3|V W | ~ 5%), could account 
for the solar neutrino puzzle. 

Solution 3 corresponds to choosing r Y = 0.4. m s /m d = 24.6 is within the allowed 
range. However, the mass ratio v T /v^ is - 3.6 x 10 3 , and sin6» MT ~ 3|V C &| so v e - 
oscillation cannot be responsible for solar MSW. As in other cases, v e - v T MSW 
oscillation with a 6% mixing is a viable possibility. 

Let us now re-instate the CP-violating phases a and (3 in the vev's perturba- 
tively. Small values of the phases are sufficient to account for realistic CP violation 
in the quark sector. We shall present details for the case of Solution 2 only, others 
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are similar. We also tried to fit all the charged fermion masses and mixing angles 
for large phases, but found no consistent solution. 

First we make a basis transformation to go from the basis where M u is diag- 
onal to one where the matrix hn u is diagonal. It is easier to introduce phases in that 
basis. For a = 3.5°, (3 = 4.5°, the CP-violating parameter J for the quark system[24] 
is J ~ 1 x 1CT 5 , which is sufficient to accommodate e in the neutral K system. The 
leptonic CP violating phases are correspondingly small, for eg., the analog of J is 
Ji ~ 7 x 1(T 5 . These small phases modify the first family masses slightly, but the ef- 
fect is less than 10%. Our predictions for the neutrino mixing angles are essentially 
unaltered. 

Recently, three more solutions for neutrino masses have been found by Lavoura, 
which give good fits to the charged fermion masses [25]. We refer to ref.25 for details 
on them. 

5. Supersymmetric minimal SO(10) and its implications: 

In this section, we consider the supersymmetric version of our model. The 
fermions are part of the chiral superfields transforming as the 16-dimensional spinor 
representation of the 50(10) group and are denoted by ip. We will choose the Higgs 
multiplets for the supersymmetric version the same way for the non-SUSY case, 
i.e., Higgs fields belonging to 10, 126 and 210 representations (denoted by H, A 
and $ respectively) but we have to include a 126 dimensional multiplet (denoted by 
A) to maintain supersymmetry down to the electro-weak scale. The most general 
superpotential for the model can be written as: 

W = W m + Mi* 2 + Ai* 3 + A 2 AA$ + M2 AA + A 3 $AH + A 4 $AH + ^H 2 (17) 

where 

W m = hablpalpbH + f ab ip a ip b A . (18) 

The gauge and supersymmetry breaking is determined by writing down the F- 
terms and setting them to zero at the GUT scale. The detailed analysis of this is 
given in ref. 26. Here, we simply quote the result of that paper which says that 
if we do not fine-tune any of the parameters of the superpotential in the above 
equation, the vanishing of F-terms requires that the 5O(10) symmetry break down 
to the standard model in one step. This implies that v R ~ M v . This has several 
interesting implications as we will see later. Unfortunately, as in all SUSY GUT 
models, in the supersymmetric limit, there are several degenate vacua corresponding 
to the symmetry groups i) 5C/(5), ii) SU(2) L x SU(2) R x U(1) B -l x SU(3) c ; iii) SU(2) L x 
U(l) ha x U{l) B -L x SU{3) C iv) SU(5) x U(l) v) SO(10). Once supergravity corrections 
are included this degeneracy disappears, and by appropriate choice of parameters, 
we can always choose the global minimum to correspond to the standard model. 

Doublet- Triplet Splitting: 
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The second point to notice is that in the Higgs potential derived from this 
superpotential, there are terms of the form AAAH proportional to A 2 A 3 and A 2 A 4 so 
that a vev for the (2, 2, 15) sub-multiplet of 126 is induced once 50(10) and B - L are 
broken, in analogy with the non-SUSY case. The detailed analysis of this question 
is tied to the question of doublet-triplet splitting of SUSY-GUT models. Let us 
therefore briefly discuss this question below. In the minimal 50(10) model there are 
eight light Higgs doublets, four H u type (i.e. capable of giving masses to up quarks 
) and four H d type which can give mass to down type quarks. After GUT symmetry 
breaking , one therefore has a 4x4 matrix and for the theory to reproduce low 
energy physics, this matrix must have two massless eigen-modes. In general this 
will require some fine tuning of parameters so that all doublets do not acquire GUT 
scale mass. At the same time, this model has a 5x5 mass matrix involving the 
color triplet fields. Successful doublet-triplet splitting requires that the same fine- 
tuning that leaves a pair (H u ,H d ) massless should not at the same time leave any 
triplets below the GUT scale, otherwise, not only will there be problem with proton 
life-time, there will also be problem with unification of couplings. The doublet- 
triplet splitting is called natural when the required fine tuning is guaranteed by 
some symmetry. 

To study this question in detail, let us isolate the various doublet and triplet 
fields in the model : The four up Higgs doublets are denoted by: H u from 10, a u 
and a u from 126 and 126 respectively and a u from 210 and similarly four down 
type doublets from the same Higgs representations. The five color triplets fields 
belonging to 3 rep. and electric charge -1/3 are denoted by Ci from 10, 7 X from 126, 
72 from 126, u from (1,3,15) in 210 and 8 R from (1,3,10) in 126. There is of course 
a similar set belonging to 3 to be denoted by (2,72? 71 > <*> an d Sr. The 4x4 Dirac 
matrix involving the doublets can be written as follows: ( the basis for this matrix 
is (H u ,a u ,a u ,a u ) denoting the columns and corresponding fields with subscripts d 
denoting the rows). 



M 



D 



f M3 




\4v1j 


A4VR 


\ 







M2 








M2 















A 2 ^K 


Ai 


/ 



(19) 



In the above matrix, v R and v v stand for the vev's of the 126 and 210 dim. 
Higgs multiplets respectively and the symbols fa denote some combination of in 
with the vev contributions. The color triplet matrix in the basis (6,71,72, for 
columns and corresponding 3 denoting rows is: 
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A 3 W(7 


A4UR 


A 4 W(7 


\ 


XiVu 















A 3 «I7 







A 2 «K 


A 2 w;7 











Mi 


A 2 «i? 






A 2 V[/ 





A 2 «i? 


M2 


/ 



(20) 



We wish to point out that, the entries in the above matrices do not con- 
tain the detailed numerical coefficients nor the various vevs corresponding to the 
multiplets- rather just the orders of magnitudes. It is clear that, one needs fine 
tuning to get a light pair of Higgs doublets. For our scheme to work, we need a 
linear combination of H u and a u (and similarly with subscript d) to remain massless 
at the GUT scale. There are several ways to achieve this. One particularly simple 
way is to set A 3 = A 4 = and add to the theory a Planck scale induced term of 
the form 10 126 126 10/ Mpi such that it balances against the ^3 term after GUT 
scale symmetry breaking. At this stage, there are two light doublets, which arise 
only from the 10 Higgs. But if we add a further Planck induced term of the form 
10 126 126 126/Mpi, then this causes an admixture of a u with H u and similarly for 
H d as required. The corresponding color triplet will however be somewhat light 
and considerations of proton decay then suggests that the coefficients of the Planck 
induced term be rather large. 

The pattern of the predictions for the neutrino masses and mixings in the 
SUSY version is very similar to the non-SUSY version. 

6. Summary and Conclusion: 

In summary, we feel that if the present indications for neutrino masses in the 
micro-, milli-, and eV range continues to receive full confirmation, then the SUSY 
50(10) may provide the correct framework for the next stage in the search for further 
unification of forces and matter. Our recent work described in this report would 
suggest that the relevant 5O(10) theory could very well be the minimal one with 
just three Higgs multiplets needed to break the gauge symmetries appropriately. 
Just like the minimal SU(5) had the unique feature that it predicted a value for 
both proton lifetime and sin 2 9 w , and could therefore be tested experimentally, the 
minimal 5O(10) grand unified model leads to definite predictions for light neutrino 
masses and mixing angles in terms of observables in the charged fermion sector 
and could therefore be tested in the next round of neutrino oscillation experiments. 
This in our opinion elevates minimal 5O(10) theory above other competing GUT 
scenerios that can lead to nonvanishing neutrino masses such as SU(5) x U(l) or E 6 
theories in terms of testability. In a sense, the status of minimal 5O(10) is better 
than the minimal SU(5) since here we fit all charged fermion masses without invoking 
additional Higgs multiplets. Of course if we want further predictability , say in the 
charged fermion sector, we will certainly need new symmetries. This line of research 
is currently being pursued. 
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We also wish to point out that our approach here has been orthogonal to 
some other recent attempts [27] based on grand unification; we have kept the Higgs 
sector as simple as possible and followed its consequences. In particular, our results 
are manifestly stable under radiative corrections. Moreover, we have not used any 
sort of family symmetry to make the model predictive. Predictivity in the neutrino 
sector follows from minimality. To repeat again, we have found three different types 
of solutions for the neutrino spectrum. In Solution 1, the v e - mixing angle is 
near the Cabibbo angle, while Solutions 2 and 3 have it much smaller. In all cases, 
v e - v T mixing angle is predicted to be near 3\V t d\ - 0.05 and - v T mixing angle is 
~ 3| Vct,| ~ 0.15 with the mass ratio m VT /m Vti > 10 3 . None of the solutions admit the 
vacuum oscillation solution to the solar neutrino puzzle. If it is due to small angle 
non-adiabatic MSW, as in Solution 2, z^-tv oscillation should be observable in the 
forthcoming experiments. 

On the theoretical side, there is one important lesson for 50(10) model builders, 
i.e. whenever, there is a 10 contributing to charged fermion masses and 126 con- 
tributing to the right-handed neutrino masses, one must include the induced vev 
contribution also in the charged fermion sector. This will certainly effect the predic- 
tivity of the models in the charged fermion sector. Our model also has the potential 
to solve the strong CP problem: note for instance that if A 4 is set to zero, we have 
a model with a softly broken PQ symmetry. 

The work of R.N.M. was supported by a grant from the National Science 
Foundation . I would like to thank K.S.Babu for collaboration in this work and Lee 
Dai Gyu for discussions. 
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